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ABSTRACT 

Hydraulic ram is the dynamic loading of fuel tanks when impacted by 
bullets or other projectiles o During impact and penetration of the fuel cell, 
intense pressure waves are generated by the projectile. 

A ballistic range was built and experimental testing was conducted 
to study hydraulic ram phenomena. A Q.22 caliber rifle was used to accel- 
erate projectiles at velocities in the range of 0.38 km/sec into a trans- 
parent, water filled tank. Shape and intensity of the shock wave pressure 
pulse induced as a result of projectile impact were determined using a dual 

shadowgraph system. Peak pressures were found to be as high as 4.50 x 
3 2 

10 kg/cm . The rate of energy transfer to the fluid by the projectile was 
determined experimentally and compared with analytical predictions. A 
characteristic time was defined, establishing the separation point between 
the shock and cavity phases of hydraulic ram. 
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I , INTRODUCTION 



When enclosed fluid containers are impacted by high speed pro- 
jectiles, severe damage can occur as a result of the phenomenon known 
as hydraulic ram. Numerous experimental investigations of the hydraulic 
ram effect have been conducted including research by Cole, Stepka , 
Lundstrom, and Yurkovich (Refs. 1-7). An investigation, presently being 
conducted at the Naval Postgraduate School, graphically demonstrates 
this effect. The damage shown in Figs. 1 and 2 occurred when a plexi- 
glas test tank, simulating an aircraft fuel cell, was impacted by a 2.92 
gram, 0.222 caliber projectile traveling at a velocity of 0.889 kilometers 
per second (2910 ft/sec). 

In determining aircraft survivability the destruction due to hydraulic 
ram becomes of paramount importance and, in fact, is the major cause of 
damage when an aircraft is hit by small arms fire or warhead fragments 
(Ref. 6). This is further substantiated by the fact that fuel tank (or cell) 
components comprise the greatest portion of the target area in modern 
attack and fighter aircraft. There are numerous examples of catastrophic 
failure of fuel cells which were subjected to ballistic impact (Ref. 6). 
Other kill modes due to hydraulic ram include: fuel starvation, explosion 
and fire, and damage to aircraft components adjacent to the fuel cell. 

In space flight the hazard of meteoroid impact into liquid propellant tanks 
is another area of concern because an impact of sufficient energy could 
cause catastrophic bursting of the tank wall. 
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Fig. 1 Hydraulic Ram Damage to the Test Tank. 



11 




12 



Fig. 2 Hydraulic Ram Damage. 




Concern with the aircraft fuel tank vulnerability problem did not 



begin with the recent air conflict in Southeast Asia. Since World War II, 
there have been efforts to reduce hydraulic ram damage to an acceptable 
level; however, with few exceptions, the efforts have been "trial and 
error" procedures. The result has been that many schemes have been 
tried, but none has been found fully effective. While not being satis- 
factory, many of these ideas have shown promise in defeating at least 
some of the hydraulic ram components. Several important fuel tank defense 
schemes were developed and tested by North American Rockwell during the 
late 1960's. These concepts are published as a survivable fuel tank sys- 
tems handbook (Ref. 8) in which the basic design features are divided into 
"intrinsic" and "parasitic" defense concepts. Intrinsic defense concepts 
utilize specific aircraft configuration considerations to reduce system 
vulnerability. Examples would be: tank shapes with smooth contours, 
low surface-to-volume ratios, few corners and fittings, and fuel system 
considerations to isolate the cell. This type of defense concept is thus 
"intrinsic" to the basic airframe design and can be effectively applied 
only when the aircraft is in preliminary design stages . Parasitic defense 
concepts include protective devices and systems added to the aircraft 
primarily to enhance fuel cell survivability; these concepts include ele- 
ments not otherwise essential to normal aircraft function. Examples are: 
self-sealing tanks and the addition of reticulated foam inside fuel cells 
to reduce the hydraulic ram effect. 
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The purpose of this study is the investigation of the hydraulic ram 



phenomenon itself. Fuel tank response to projectile penetration cannot 
be fully understood without a prior comprehension of the individual 
hydraulic ram components. This understanding, plus a knowledge of 
events which occur during the damaging process, must be obtained in 
order to change the tank arrangement in such a way as to properly influ- 
ence the damage. Once fuel cells are protected against this type of 
threat, overall aircraft vulnerability will be lessened. 
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II. NATURE OF THE PROBLEM 



It has been found (Refs. 6 and 7) that the loading and deformation 
of fuel tanks by hydraulic ram is a combination of a number of events . 
These events, known as the hydraulic ram components, are outlined in 
Figs. 3 and 4. The column on the left in Fig. 3 describes the motion of 
the projectile as it passes through the tank. When the projectile enters 
the tank, a shock wave is formed and cracks may develop due to the high 
pressure and stresses near the entry point. As the projectile traverses 
the fluid, high pressure is generated, and energy is imparted to the fluid 
through projectile drag. This increased energy sets the fluid in motion 
forming a cavity. Both the fluid motion and the subsequent cavity forma- 
tion and collapse impose stresses on the tank walls. Finally, if the pro- 
jectile has sufficient kinetic energy, it v\?ill exit the tank. As it exits 
from the tank, a local area of compressed fluid is formed, and cracks are 
produced around the exit hole. As shown in Fig. 4, the wall loadings 
plus the cracks at the entry and exit points can result in various types of 
damage . 

Excluding crack formation due to the actual entry and exit of the 
projectile, hydraulic ram is composed of two main separable effects; a 
shock phase and a drag or cavity phase, both of which are important in 
causing damage. During the initial impact, the projectile impulsively 
accelerates the fluid, thus generating an intense pressure field of short 
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Fig. 3 Outline of Hydraulic Ram Phenomenon 
adapted from Reference 10. 
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Fig. 4 Hydraulic Ram Damage Sequence 
adapted from Reference 10. 



duration. This phenomenon is characteristic of the shock phase of 
hydraulic ram. The contribution of the shock-phase pressure pulse 
toward fuel cell failure is particularly important in the case of small, 
high velocity projectile impact. After passage of the projectile, pressure 
phenomena become increasingly complex due to: 1) multiple sources of 

the pressure, 2) reflections of the pressure waves, 3) motion of the 
fluid, and 4) tank wall deflection. These events lead to the formation 
of a cavity and mark the beginning of the cavity phase. The formation of 
a large cavity behind the projectile is the most striking phenomenon 
associated with ballistic penetration of fluids. The pressure field pro- 
duced by the cavity formation and collapse is of low amplitude but long 
dux'ation and has a dynamic behavior loosely analogous to that of a cavity 
formed by an underwater explosion. 

The shock and cavity phases that have just been described have 
been treated in this analysis as two separate effects. A clearer under- 
standing of the relationship between the shock and cavity phenomena can 
be gained by considering the projectile energy transfer time-history as 
shown in Fig. 5 (Ref. 9). Two separate regions can be distinguished 
which correspond to the two main effects present in hydraulic ram. The 
first region is characterized by a high energy transfer from the projectile 
to the fluid. It is while the projectile is in this region that the shock 
wave is formed and begins to propagate. The second region is character- 
ized by a lower and nearly constant rate of energy transfer to the fluid . 
This second region, which lasts for the remaining time that the 
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ENERGY TRANSFERRED TO THE FLUiD 




TIME 



Fig. 5 Energy Transferred to the Fluid by a 
Projectile Versus Time. 
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projectile is traversing the tank, is known as the drag or cavity phase. 

In this phase, projectile energy is transformed into the kinetic energy of 
fluid motion which leads to the formation of a cavity. 

The work here was done to further the development of analytical 
methods used in understanding the complex hydraulic ram phenomenon. 
Once an ability to predict the response of fuel cells to ballistic impact 
exists, overall aircraft survivability will be enhanced. 

A. SHOCK PHASE 

Shock wave formation and propagation due to projectile penetration 
into a fluid-filled tank will result in very high local pressures which can 
cause catastrophic failure of the entry wall. This result has been shown 
experimentally (Refs . 2 - 5) , and it provides the motivation behind the 
study of the shock phase of the hydraulic ram phenomenon. 

The shock or impact phase can itself be regarded as a succession 
of interrelated events which will be discussed in their order of occurrence. 
When a fuel tank wall is impacted by a high velocity projectile, dynamic 
stresses are generated and transmitted through the wall. These stresses 
are due to the puncturing and cratering action from the projectile impact 
and take the form of compressive and circumferential bending stress in 
the wall (Ref. 3). As the projectile is penetrating the tank wall, it 
impulsively accelerates the fluid directly behind the impact point. An 
instant before impact the fluid is at rest, and in the next instant it is 
traveling at the velocity of the projectile. This sudden acceleration of the 
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fluid by the projectile during impact and penetration generates an intense 



pressure field bounded by a shock wave. Tank wall stresses are induced 
by this pressure field which, when added to the dynamic stresses due to 
cratering, can cause catastrophic failure of the entry wall. The shock 
wave pressure field expands in a hemispherical shape away from the 
impact point with a velocity greater than the acoustic speed of the fluid. 
Both the velocity of the shock front and the intensity of the related pres- 
sure field vary directly with the amount of energy transferred to the fluid 
by the projectile during the initial impact (Refs. 3 and 4). The amount 
of impact energy which is transferred to the fluid is dependent on a num- 
ber of parameters including; projectile size, density, and geometric 
shape; tank wall material and construction; and the physical properties 
of the fluid (Ref. 5). In general, for any given time after impact, small 
high-velocity projectiles lose a larger part of their impacting kinetic 
energy than more massive, low velocity projectiles (Ref. 2). This 
observation indicates that for the same impact kinetic energies, a more 
damaging pressure pulse for wall fracture will resirlt from impact by 
small high velocity projectiles than with more massive low velocity 
projectiles . 

Although the peak pressure due to projectile impact is high, the 
pressure field is rapidly weakened because of the geometric expansion 
of the shock wave. An experimental example of this pressure attenuation 
is provided by Stepka (Ref. 2). Impact with a 0.9 gram steel cylinder at 
an impact velocity of 4.27 km/sec resulted in pressures of 4.63 
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G 2 

giganewtons per square meter (0.67 x 10 Ib/in ) at a distance of 1 . 9 
centimeters (3/4 inch) from the impact point. However, the pressure at 
the shock front decayed to ambient values within 13 centimeters (5 
inches) from the impact point. These results indicate that even for im- 
pacts with high velocity and high kinetic energy projectiles, the side or 
rear walls of any large tank would not be subjected to significant pres- 
sures due to the shock phase. 

In summary, the shock phase of hydraulic ram is characterized by 
an intense pressure field of short duration. It is essentially a "point 
source" energy release and results in an expanding, hemispherically 
shaped, pressure wave front that emanates from the point of impact. The 
kinetic energy of the projectile is the most important parameter in deter- 
mining the speed of propagation and strength of the shock front. The 
shock-phase pressure pulse contributes to fuel cell failure mainly in the 
case of small high velocity projectile impact. 

B. CAVITY PHASE. 

After the initial impact with the tank, which is responsible for the 
shock phase, the projectile continues to decelerate as it traverses 
through the fluid . It is during this time that the cavity phase of the 
hydraulic ram phenomenon is generated . 

When an object moves through a fluid, it feels a resistance known 
as "drag.." Drag forces are due to an uneven pressure distribution and 
viscous effects of the fluid. If the pressure distribution is integrated in 
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the direction of motion, the resulting force is known as "pressure drag." 
Compared with the drag due to viscous effects, pressure drag is the pre- 
dominant drag force (90%) at projectile velocities (Ref. 10). 

Drag is usually expressed in terms of a nondimensional "drag 
coefficient," which is dependent on the density of the fluid, as well as 
on the size, velocity, and geometric shape of the object. A typical drag 
coefficient for a high-velocity projectile is 0.15, but this value will 
increase by a factor of ten or more for a tumbled or flattened projectile 
(Ref. 10). In the analytical study of hydraulic ram, this coefficient can 
be used to account for the much greater pressure generating capability of 
deformed and tumbled projectiles. 

Applied to the projectile surface, the drag forces generate a pres- 
sure gradient between the projectile and the fluid. This causes the sur- 
rounding fluid to move away from the surface of the projectile. In this 
manner, the projectile energy is gradually transformed into the kinetic 
energy of fluid motion. In contrast to the shock phase, the fluid is 
accelerated gradually rather than impulsively so the peak pressure is 
lower. However, the duration of the pressure pulse is much longer. 

Once a sufficient amount of energy is transferred to the fluid, it will 
form a void, known as a cavity, around and behind the projectile. The 
cavity which is formed would be a complete vacuum except for two 
phenomena. It is filled with fluid vapors from pressure vaporization 
effects and with entrained air from the projectile entry and exit points. 
After passage of the projectile, the hydraulic ram cavity continues to 
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expand until all fluid momentum has been dissipated. The momentum is 



lost to the surrounding fluid by fluid pressures which resist cavity growth. 
These same pressures including pressure due to the force of gravity 
cause the cavity to collapse. Collapse of the over-expanded cavity 
would be exactly opposite to cavity growth except for the presence of 
fluid vapors and trapped air. These vapors and air act as a spring to make 
complete collapse impossible. From this point on the cavity acts much 
like an explosion cavity; the cavity oscillations generate low amplitude 
but long duration pressure pulses. The cavity may oscillate through 
several cycles until the energy of the fluid is dissipated. The pulse width 
and peak pressure produced by the cavity oscillation depend on the total 
energy transferred to the fluid, the amount of gas trapped in the cavity, 
and the shape of the cavity itself (Ref. 6). Pressure pulses due to the 
cavity phase of hydraulic ram are generally more damaging than those of 
the shock phase due to their longer duration. 

The rate of energy transfer to the fluid is the major factor that 
determines the severity of hydraulic ram. Projectile kinetic energy dis- 
sipation is proportional to the square of the velocity decay as the pro- 
jectile traverses the fluid. This indicates that the projectile deceleration 
time-history as it passes through a fluid is of great importance in the 
study of hydraulic ram. 
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III. .ANALYTICAL APPROACH 



TO CAVITY OSCILLATION AND PRESSURE FIELD 



Two analytical problems were investigated concerning the oscilla- 
tory cavity phase of hydraulic ram. The first problem involved the deriva- 
tion of a method for predicting the collapse time for a cavity when the 
vapor pressure within was negligible. The assumption applied to this 
problem was that the cavity contained a vacuum. The second problem 
can'ied the same basic analysis even further by modeling the oscillation 
of a cavity with finite internal pressures. This cavity was found to behave 
much like a spring in its oscillatory pattern. Isentropic compression and 
expansion were the basic assumptions used in the second problem. 

A. COLLAPSE OF A CAVITY CONTAINING A VACUUM VERSUS TIME 

Potential flow relationships for an incompressible fluid were used 
in the formulation of this problem. The cavity vs/as assumed to be spheri- 
cal in shape with an initial radius Rq. R = R(t) describes the motion of 
the fluid -gas surface which travels at the velocity u. From continuity, 
the flow field is described by 



where r is the instantaneous cavity radius . Potential flow theory 
along with the momentum equation for an incompressible fluid was 




( 1 ) 
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used to obtain the equation of motion for a cavity. It v^as found 
(Appendix A) that; 




r 



( 2 ) 



where P (R, t) was the pressure in the fluid at the cavity wall. Assuming 
that P (R, t) was Icnown, equation (2) was manipulated into the general 
form 




For the first problem with R = Rq , Rq = 0 and P (R, 0) = 0 , the 
integral yields the time for the cavity to reach a fraction of its initial 
radius 




/T-x' 



(4) 



where X = R / R^ . This equation was numerically integrated by compu- 
ter. The result describes the motion of a spherical cavity v.'ith time 
when the cavity contains a vacuum and starts at some initial radius Rq • 
Solutions for two cavities with varying initial conditions are shown in 
the computer outputs of Figs. 6 and 7. These plots give valuable insight 
into the frequency response and an order of magnitude for the time to 
cavity collapse. Figure 6 gives a time to collapse of 0.011 seconds 
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R/R 




Time (seconds) 

Fig . 6 Cavity Collapse Time from an 
Initial Radius of 0.5 Foot. 
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R/R, 




Time (seconds) 

Fig, 7 Cavity Collapse Time from an 
Initial Radius of 1.0 Foot. 
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2 

for a cavity with the parameters: Rq= 0.5 ft, Poo = 2200 Ib/ft , and 

3 

y^= 1.94 slug/ft , P’igure 7 is the solution to the problem with the 

2 

following initial conditions: Rq ~ ^*0 ^oo = 2500 Ib/ft , and 

1.60 slug/ft^ . 

B. ISENTROPIC CAVITY OSCILLATIONS 

The second problem, which assumed the cavity was filled with gas 
and followed a series of isentropic compressions and expansions, was 
approached in a manner similar to the first problem. Starting v/ith 
R ( 0) = Rq 2 nd P (R, 0) = Pq , it was known that for isentropic 
compression 




which when applied to a homogeneous gas-filled sphere yields 




This relationship was derived by assuming that R << a, where a was the 
speed of sound in the gas. Combining equations (3) and (6) it was found 
that 
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A dimensionless time was defined as 




/ 




( 8 ) 



and combined with equation (7). The final expression for cavity motion 



v>/here c< = Pq/( V ~1) Poo and X = Rq/R. Equation (9) was solved 
through the use of a numerical integration computer program. The cavity 
radius, represented by X = Rq/R, was found to oscillate about the equi- 
librium position X = 1. This response was expected, since the cavity 
started at an initial radius Rq , and there was no damping term included 
in the analysis. This departure from the true cavity response was brought 
about by the assumption of isentropic compression and expansion of the 
air inside the cavity and the neglect of compressibility effects in the 
fluid outside the cavity. However, valuable information can be gained 
concerning the frequency and peak pressures associated with the cavity's 
repeated expansion and collapse. To aid in the study of these parameters, 
the computer output was designed to reflect the cavity pressure fluctua- 
tions with time. This was accomplished by plotting the output in the 
following form: cavity pressure/ambient fluid pressure versus time. This 




(9) 
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pressure is directly related to the cavity radius as shown by equation 
(6). An example of the output, which was computed for an initial pres- 
sure ratio of 6.67, is shown in Fig. 8. This plot illustrates the cavity 
pressure oscillation with time, and further it predicts that the peak 
pressure occurs when the cavity is at its minimum radius. This fact 
has been shown experimentally by Cole (Ref. 1). 



31 



p/p, 




t ■ 



Fig, 8 Hydraulic Ram Cavity Pressure Oscillation. 
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IV. APPARATUS 



The test apparatus consisted of projectile accelerators, test speci- 
mens, a test tank, and associated instrumentation for evaluating the 
effect of high-velocity impact into liquid filled containers. 

A. PROJECTILE ACCELERATORS AND PROJECTILES 

A 0.22 caliber rifle was used to accelerate the projectiles. The 
rifle was mounted on a stand, as shown in Fig. 9, and located about 6.9 
meters from the target tank. Impact tests were made with projectile 

I 

velocities in the order of 0.381 km/sec (1250 ft/sec). 

The cartridges were factory loaded 0.22 caliber long rifle ammuni- 
tion with a solid lead projectile weighing 1.95 grams. 

B. TEST TANK 

The test tank used for the impact tests was fabricated from trans- 
parent Plexiglas, Type G, 25.4 millimeters (1 inch) thick. The tank was 
square in shape, with the plexiglas making up the top, bottom, and two 
sides. The front and rear of the tank had removable and replaceable 
walls into which impact was made. A photograph of the tank with its 
support table is shown in Fig. 10. In this investigation, the removable 
walls were made from 7075-T6 aluminum sheet 1.27 millimeter (0.050) 
thick, secured to the transparent walls by a clamping frame. 



33 



Jf- ' 




34 



Fig. 9 Ballistic Range Rifle Mount 




Fig. 10 Test Tank Installation . 
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C. RANGE LAYOUT 

The basic elements of the ballistic range are shov\n in schematic 
form in Fig, 11. The down-range photograph (Fig. 12) shows the velocity 
sensors and the three shadowgraph stations located along the flight path. 
These shadowgraph stations were used in conjunction with the velocity 
sensors to accurately determine the velocity and attitude of the projectile 
passing down the range, 

D. INSTRUMENTATION 

The projectile was a chronograph screen with a 5-volt d.c. potential 
across it. When the screen was broken by a projectile, a pulse was sent 
through a signal conditioner, then to the delay unit, and finally to a 
counter. Two Monsanto 10 IB counters were used for velocity measure- 
ment. The delay units were used to trigger the shadowgraph spark 
sources at the proper time. Figure 13 shows the circuitry of the two 
delay units used in timing the dual shadowgraph system of the test tank. 
The major component of each delay unit was a Signetics N74123 retrigger- 
able, monostable multivibrator. These units generate variable time delays 
from 200>cs to 1700 a<s and from to 18^s / respectively. A spark 
source triggering capacitor was built into each delay unit. The signal 
conditioner, shown in Fig. 13, was needed to compensate for the wide 
range of resistances found in the chronograph screens. The signal 
conditioner was basically a transistor switch using a 2N3404 transistor. 
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Fig. 11 Ballistic Range Layout, 
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Fig. 12 Ballistic Range Configuration (Down- range View) 



I 




Vcc ( 5 Volts d.c. ) 




Fig. 13 Delay Unit Circuitry. 
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E. SHADOWGRAPI-I SYSTEM 

The shadowgraph system associated with the test tank consisted 
of two similar stations . Figure 14 is a schematic representation of the 
apparatus showing a down-range view. Each shadowgraph station was 
composed of a spark source, a 16-inch collimating mirror, and a 4-x 5- 
inch photographic plate. The spark source was located at the focal point 
of the mirror to provide the necessary collimated light. Polaroid Type 57 
high-speed film was used to record the shadowgraphs through the water 
filled test tank. The complete shadowgraph system is shown in Figs. 15 
and 16. 
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station 1 




Fig, 14 Schematic of Shadowgraph Apparatus 
(Down- range View) . 
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Fig. 16 Dual Shadowgraph System 

(Plane Mirror and Film Holders). 
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V. RESULTS AND DISCUSSION 



This research was conducted as a preliminary investigation of the 
mechanisms involved in hydraulic ram. Primarily, the phenomena studied 
were the characteristics of the shock wave pressure pulse generated in 
water from impact by a projectile and the rate of energy loss of the pro- 
jectile in the water. By the characteristics of the shock wave pressure 
pulse one means the shape and progress of the shock wave with time 
and the magnitude of the pressures at the wave front. 

A. SHAPE OF THE SHOCK WAVE FRONT 

The investigation of Ref. 3 indicated that the shape of the pressure 
wave front generated in the water by a projectile impact was hemispheri- 
cal. However, the data of this reference were primarily limited to impacts 
with spherical projectiles. Since the typical threat to an aircraft fuel 
tank has the standard bullet shape, the projectiles used in the investiga- 
tion were solid lead bullets . The influence of lead bullets , which deform 
as they enter the tank and traverse the fluid, was explored and will be 
discussed in a later section of this report. 

The shape and propagation of the pressure wave front resulting from 
the projectile impact were obtained from the dual high-speed shadowgraphs. 
A typical dual shadowgraph sequence of the pressure wave front progress 
is shown in Fig. 17. The figure illustrates the results of the impact by a 
1.95 gram lead bullet at a velocity of 0.374 kilometers per second 
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(1226 ft/sec) into water in the transparent plexiglas test tank. For this 
test the projectile passed through a 0.5-inch diameter, "pre-punched" 
hole in the entry wall. A strip of masking tape was used to cover the 
hole and contain the water. A hemispherical wave front was produced, 
and the center of the hemisphere remained at the point of impact as the 
wave front expanded. The photographs of Fig. 17 show that there was 
no significant effect of the forward momentum of the projectile to move 
the wave front in the direction of flight. 

A commonly used model for analytical prediction of the hemi- 
spherical wave strength and motion is the blast wave theory. Blast 
wave theory assumes sudden release of energy centered at the projectile 
penetration point. The behavior of the fluid properties as a function of 
radius will not be consistent with that model. See, for example. 

Chapter 4 of Sedov (Ref. 11). 

An alternate and more accurate model would be to assume an 
expanding hemispherical piston due to projectile penetration. The rate 
of change of piston radius would be 




( 10 ) 



The quantity ()A// () / equals the rate at which the projectile dis- 
places fluid volume during penetration. 
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27 after impact 




32 after impact 



Fig. 17 Shock Waves Resulting from a 
Pre-punched Wall Impact. 
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Tests were also performed with a solid front wall on the test tank. 



The entry walls used were 7075 -T6 aluminum sheet 1.27 millimeter 
(0.050 in.) thick. An example of the resulting wave front due to solid 
wall penetration is shown in Fig. 18. This set of photographs illustrates 
the obscuration of the interior of the shock wave due to the projectile 
penetration of the entry wall. The shadowgraphs show that the hemispher- 
ical wave front was produced as in the "pre-punched wall case." The 
other wave front in the water, which appears as a straight line oblique 
to the wall, was caused by the propagation of shear and dilatation waves 
(Ref. 12) moving through the metal at a faster rate than the waves through 
the water. The rate of propagation of this oblique wave was found to be 
proportional to the bulk sonic velocity (dilatation wave) for aluminum 
(6 . 2 km/sec) . 

B. PROGRESS OF THE SHOCK WAVE WITH TIME 

Impacts with projectiles of various kinetic energies were made to 
determine the influence of kinetic energy on shock wave velocity. The 
average velocities of the shock waves were determined from shadowgraphs 
similar to those in Figs. 17 and 18. The shadowgraphs were used to ob- 
tain time and position information for the shock wave as it expanded and 
propagated through the fluid. The slopes of the resulting time - distance 
plots supplied the average propagation velocities of the shock waves. 

Figure 19 is a plot of the resulting shock front velocities for an 
impact kinetic energy of 135.2 joules. This plot shows that the shock 
front velocity decays to the acoustic speed in water by a distance 
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18 y£cs after impact 




25 yus after impact 



Fig. 18 Shock Waves Resulting from a 
Solid Wall Impact. 
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SHOCK VELOCITY (km/sec) 




Fig. 19 Shock Front Velocity Generated in Water 
by Projectile Impact. 
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of about 11 centimeters (4.3 in.) for this impact energy. The acoustic 
speed for water used in this investigation was 1.5 kilometers per second 
(4920 ft/sec) . 

C. PRESSURES GENERATED IN THE FLUID 

The pressure generated in the fluid by the projectile was divided 
into two regions. The first region consisted of the pressure associated 
with the shock front propagation and was considered as the shock phase 
of the hydraulic ram phenomenon. The second pressure region was con- 
centrated in the immediate vicinity of the projectile itself. 

1 . Shock Front Pressure 

Pressures generated in the water at the shock front correspond 
to the measured shock front velocities and were determined from data given 
in Table 2 . 2 in Ref. 1. Figure 20 shows the pressures generated at the 
shock front corresponding to the data given in Fig. 19 of this report. The 
result shows that extremely high pressures are generated in the water 
near the impact point. For example, impact with a 1.9 5 gram bullet at an 
impact velocity of 0.372 kilometers per second (1222 ft/sec) resulted in 
a pressure of 3750 kilograms per square centimeter ( 53,335 Ib/sq 
in.) at a distance of 2.24 centimeters (0.88 in.) from the impact point. 
These pressures decay rapidly due to geometric expansion of the shock 
front and, as in the case of the shock velocity, reach ambient values in 
a short distance . 
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Fig. 20 Shock Front Pressure Generated in 
Water by Projectile Impact. 
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2 . 



Pressure Field Around the Projectile 



As with any object moving through a fluid, the projectiles 

created a dynamic pressure in the water due to the stagnation conditions 

on their front surface. Dynamic pressure is related to the density of the 

fluid and the square of projectile velocity. For the velocities of 0.372 

kilometers per second used in this research, pressures of 710 kilograms 

4 2 

per square centimeter (1.01 x 10 Ib/in ) were found to exist in the vicin- 
ity of the projectile. 

Further insight into the pressure field around the projectile was 
gained by examining the light diffraction patterns illustrated in Fig. 21. 

The concentration of light at the front surface of the projectile was caused 
by a change in the refractive index of the fluid due to the large increase 
of pressure and density in this region. Light is refracted toward the larger 
refractive index which accompanies higher density and pressure in the fluid 
A qualitative analysis of the entire pressure field around the projectile can 
bo made by studying the amount and direction of deflection of the grid lines 
in this area . 

Figure 21 also shows the air filled conical cavity behind the 
projectile. This cavity is a result of the dynamic pressure produced in the 
fluid by the projectile. 

D . PROJECTILE ENERGY DECAY AFTER IMPACT 

Measurements were made of the progress of the projectile through the 
water after impact. From these measurements the projectile velocity and 
energy loss were determined. 
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Fig. 21 Light Diffraction Pattern Around the Projectile. 



One set of the shadowgraphs from which the projectile velocities 
were calculated is shown in Fig, 22. These experimentally derived veloci- 
ties and energies are plotted as a fraction of their impact magnitude 
versus distance into the tank, as shown in Figs. 23 and 24. 

A detailed analysis of the velocity decay or energy loss of a projec- 
tile traversing a fluid would be extremely difficult because of the many 
variable parameters. Assumptions must be made to account for the changes 
in the projectile's shape, drag coefficient, and material properties. Because 
of these difficulties, a simplified analysis of the projectile velocity and 
energy loss was made and compared with the experimental data. 

The following drag equation was used: 






( 11 ) 



Rearranging terms , knowing V = dx/dt, gives 



C/) /I Jx 

V z m 



(12) 



Assuming a constant deformed projectile area of 0.5 square centimeters, 
constant drag coefficient of 1,0, and constant projectile mass, integration 



yields: 



V„ 



= e 



_ FfCnAx 
zm ^ 



(13) 



Squaring both sides of this relation gives the ratio of the projectile energy 
after impact to the projectile impact energy. These calculated ratios for 
the projectile velocity and energy are compared with the experimental data 
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100 yus after impact 




264 after impact 



Fig. 22 Projectile Progress into the Tank. 
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V/ Vq= e 



2 m 



V/Vo 




Fig. 23 Projectile Velocity Decoy Ratio as a 
Function of Distance into Tank. 
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Fig, 24 Experimental and Calculated Ratios 
of Projectile Energy. 
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in Figs. 23 and 24. This comparison indicates that in spite of the sirnpli-* 



fying assumptions, the agreement between the calculated and experimental 
values of the velocity and energy after impact v/as fairly good. 

In the experimental tests, the actual projectiles underwent a con- 
tinual deformation during deceleration in the tank. An increase in frontal 
area accompanied the deformation, resulting in an average final projectile 
area of 0.5 square centimeters. This was the area used in the calculated 
approximation of the velocity and energy decay and constitutes the major 
cause for error. 

Energy decay of the projectile can be directly related to the amount 
of energy being transferred to the fluid by the projectile. More importantly, 
the slope of the projectile energy decay curve gives the rate of energy 

transfer to the fluid, assuming no projectile deformation. This energy 

\ 

transfer rate has been established as a major factor in determining the 
severity of hydraulic ram (Ref. 7). The average energy transfer rate for the 
projectiles tested was 20.6 joules per centimeter of fluid traversed. 

E. PROJECTILE DEFORMATION AFTER IMPACT IN FLUID 

The projectiles were recovered after each experimental test and were 
found to display a unique pattern of deformation. The deformation phenome- 
non was very consistent for the projectiles impacting the fluid through the 
pre-punched plates. Three examples of projectiles after this type of impact 
are shown in Fig. 25. The extent of deformation was found to vary directly 
with the impact velocity; however, all of the projectiles exhibited an iden- 
tical pattern of having one flat side and one undeformed side. This pattern 
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Fig. 25 Projectiles After Impact Tarough 
Pre-punched Entry Wall 
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of deformation suggests a positive dynamic stability for the projectiles 
as they traverse the fluid after impact. Once initially flattened, the pro- 
jectile maintains a fixed attitude due to its stability, and all further 
deformation takes place on the same flattened side. As shown in Fig. 26A, 
the dynamic stability of the projectile is maintained through an equilibrium 
condition between the pressure force, located at the center of pressure, 
and the inertia force, acting through the center of gravity of the projectile. 
Figure 26B demonstrates the effect of displacing the projectile from the equilib- 
rium position. The center of pressure moves upward on the slightly curved 
front surface of the projectile with the stagnation point; this creates a 
restoring moment which returns the projectile to the stable position. 

The projectiles which penetrated a solid entry wall before traversing 
the fluid did not demonstrate the previously discussed deformation phenome- 
non. Two typical projectiles of this type are shown in Fig. 27 . The 
small aluminum cap located above each projectile in the photograph is the 
piece of the entry wall punched out by the projectile as it passed through 
the aluminum plate. A photograph of an entry wall after a test shot is 
shown as Fig . 28 . 

F. SHOCK AND CAVITY PHASE SEPARATION POINT 

When the projectile first penetrates the wall, a strong shock wave 
pressure pulse is created. As shown in Fig. 29, this pressure decays in 
sti'ength as' it progresses outward. At time t;^ in Fig. 29 the pressure 
behind the shock greatly exceeds the stagnation pressure on the projectile 
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(a) At eq\.n librium 
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(B) Displaced from cqxxilibrium 



Fig, 26 Dynamic Stability of a Deformed 
Projectile in Water After Impact. 
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Fig. 27 Projectiles After Impact Through 
Solid Entry Wall. 
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Fig. 28 Solid Entry Wall After Impact. 
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Fig. 29 Evolution of Pressure in the Transition from Shock to Cavity Phase. 
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nose. Comparing pressures at t^ and , it is apparent that both the 
pressure behind the shock and at the projectile forvv^ard stagnation point 
are decaying; hov/ever, the decay of the pressure behind the shock is much 
more rapid. At time t^^ the stagnation pressure exceeds the pressure behind 
the shock wave. Therefore, at some point between t and tQ the two 
pressures, shock and stagnation, are equal. One could define this as 
a characteristic time separating the shock and cavity phases of hydraulic 
ram. The characteristic time for the impact energy of 135.2 joules shown 
in Fig. 29 is 32 microseconds. 

A comparison can be made to further understand the differences between 
the shock wave and projectile pressures. The stagnation pressure created 
on the nose of the projectile as it passes through the fluid can be modeled 
by steady flow and the equation (14) from Ref. 13, 







(14) 



while the discontinuity of the shock wave requires the use of the following 
unsteady flow equation (Ref. 13). 







y^i / 



Co 









y- 1 

Y^l 



y 



(15) 



For the same fluid velocity (V/ Cq = 0.5) the pressure increase is 19 per 
cent for steady flow and 84 per cent for unsteady flow. This comparison 
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explains the large initial pressures of the shock wave as well as their 



rapid decay since unsteady flow is transitory. 
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VI. ALTERNATE APPROACHES TO THE PROBLEM 



As a result of the work reported here, the following apparatus and 
test procedure are recommended for use in the future study of the hydraulic 
ram phenomenon. 

A. HIGH SPEED FRAMING CAMERA 

The use of a high speed framing camera would be the optimum way to 
visually evaluate the hydraulic ram phenomenon. A camera capable of 
recording 200,000 frames per second would be ideal for further study of 
the short-lived events which occur. This rate, a frame every 5 microseconds, 
would supply a complete time-history for the motion of the shock front and 
projectile through the fluid. 

B. PRESSURE TRANSDUCER 

Pressure transducers are essential for the complete analysis of the 
complex hydraulic ram pressure field and are recommended for use in future 
work. The pressure distribution is extremely complex due to tank wall 
reflections and the continued expansions and collapses of the cavity. 
Recording apparatus for these pressure phenomena require a high frequency 
response (about 20 kHz for cavity phase and 1 MHz for shock phase) and 
are very sensitive to location (Ref. 6). However, it is important that this 
type of direct pressure readout be compared with previously obtained 
information . 
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C. STRAIN GAUGES 

After projectile penetration, a bending moment is produced in the 
entry wall and propagates with the expanding shock v/ave pressure pulse. 
Strain gauges would be valuable in determining both the magnitude and 
propagation speed of the bending moment and its associated pressure wave 
front . 

D. HYDROGEN BUBBLE DISPLACEMENT 

In the cavity phase of hydraulic ram, knowledge of fluid motion 
within the tank is important. Hydrogen bubble flov/ visualization techniques 
are well suited for investigations requiring determination of the velocity 
field over an extended region of space in a time-dependent flow. This 
technique could be used with the transparent test tank and high speed 
framing camera to determine the pressure distribution surrounding the cavity 
as it expands and collapses. 

E. HOLOGRAPHIC INTERFEROMETRY 

In Order to study a material under shock loading, it is necessary to 
measure both the shock transit time and the profile of the transmitted stress 
wave. Wall displacement information of this type can be obtained most 
efficiently through the use of holographic interferometry. The shocked 
specimen's free surface motion measured by this technique could be studied 
and related to overall fuel tank response to hydraulic ram. 
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VII. CONCLUSION 



The following conclusions were obtained from the investigation of 
the shock wave and cavitation effects of the hydraulic ram phenomenon: 

A ballistic range and a simulated fuel tank test apparatus have been 
designed, built, and made successfully operational. 

The nonlinear equations describing cavity oscillation and collapse 
have been programmed and numerical results obtained. 

The rapid energy transfer of the shock phase produces piston- like 
energy release resulting in a hemispherically shaped pressure wove propa- 
gating from the point of projectile impact. The hemispherical shape of the 
shock wave was found to exist for impacts through thin solid metal and 
pre-punched tank entry walls . 

Very high pressures were produced in the fluid as a result of the pro- 
jectile impact. A pressure of 3750 kilograms per square centimeter 
2 

( 53,335 Ib/in ) at a distance of 2 . 2 centimeters from the impact point 

was produced by a 1.95 gram projectile with an impact velocity of 0.375 
kilometers per second. 

Good agreement was found between simplified analytical predictions 
and experimental data for energy lost by the projectile as it decelerated in 
the fluid . 

Large dynamic pressures were produced on the projectile as it passed 
through the fluid. These pressures, in the order of 700 kilograms per square 
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centimeter (10 Ib/in ), caused substantial deformation of the solid 
lead projectiles. 

A characteristic time was defined which establishes the point of 
separation between the shock and cavity phases of hydraulic ram. For an 
impact energy of 135.2 joules the characteristic time was found to be 32 
microseconds . 
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APPENDIX A 



DERIVATION OF THE EQUATION OF MOTION 
FOR THE HYDRAULIC RAM CAVITY 



From continuity one obtains 




which integrates to the form 

CL - r (i) 



Since there is potential flow 

^ ^ ^ _ R R 

^ ' J7 - ~F^ 

Integration yields 

Differentiating with respect to time gives 

j(j> _ _ R * ZRR’' j6 
^ i l~ J z' 

For an incompressible fluid the momentum eguation in terms of potential 
functions is 



^ (P P 

c)z^ V 
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Considering boundary conditions at infinity ^ Qp / c) ^ = 0 

Poo 

and C (t) = ~7T~ • Combining equations leads to 






_P _ ^ a R"' R* 

' r r 2 r"" 



For conditions inside the cavity, r = R , and the preceding equation 



becomes 



• • 

RR 




P(R.i) - Po. 



which is the equation of motion for a cavity. 
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COMPUTER PROGRAM 



//HCLPl JCB (C738, L028,AR22 ) , 'HCLP' 
// EXEC FGRTCLGP, REGION. GC=LOOK 
//FCRT .SVSIN CD * 



C/iVITY FROBLER NUREER 1 



:}f ^ a}? 5}t 5^ sit 5^ >;< 5^ ^ y: 5;^ lit >> A ^ y* 

5}t 

^ TFJS PRCBLEM USES THE NUMERICAL INTEGRA! ICN SLB- * 

•1- RCUTINE QATR TO SOLVE THE NCN-LINEAR CAVITY EGUATICN * 

* OF ROTICN. PRCbLER NUMBER 1 IS RESIGNED TC CCRPUTE 

* TFE COLLAPSE TIME CF A SPHERICAL VACUUM FILLEC HYC- 
RALLIC RAM CAVITY. 

4 :1c 

^ .* 

4 ^ 4 :4 ^ 4c c<c :{c ^ ^ :>fc ^ ^ ic Cr ^ sjt A ij; 5}; « * 5|; >S 5) ^ ^ ;f: i: :]c ^ ^ ijc :1c io{: :J; 



EXTERNAL FCT 

DIRENSICN T (200) ,XX( 2C0 ) ,AUX(20C) 

NPITE (6,.12) 

TFE INITIAL CCNDITIGNS FCR PRESSURE, DENSITY, AND 
CAVITY RADIUS ARE DEFINED. 

PINF=25C0. 

RFC = 1 .6C 
R= 1 .0 

A =- 1 . 0/ ( ( 2 . 0/ 3 . 0 ) =*■ ( P I NF / RHO ) J * C . E 

INTEGRATION PARAMETERS FOR SLEiRCUTINE QATR ARE 
ASSIGNED. 

XL=C.99CC0C00 

XU=C.99CCC000 

EFS=0.0C1 

N=lCO 

NDIR=AO 

DC 1 K=1,N 

CALL GATR(XL,XU,EPS,NCIR,FCT,Y, I EP,AUX) 

E = R’^A 
T (K )=-B^Y 
XX( K)=XL 

XL = XL-C .OICGOOOO 
1. CONTINUE 

XX( 1) = 1 .OCCO 
XX( 100) =C.COOC 
DC 2 R=l, 100 

WRITE (6,1A) M,XX(R),T(M) 

2 CCNTINUE 

SCLUTICN ARRAY IS PREPARED FOR PLOTTING USING SUE- 
RCUTINE CRAW. 



REALMS TITLE(12) 

READ (5,16) TITLE 
REAL LABLE/' '/ 

CALL DRAW ( 100 , T , XX , 0 , 0 , L A B L E , T I TL E , 0 . 0 , C . C , C , 1 , 0 , C , 5 , 
6,0, LAST ) 

C 

12 FORMAT (' I', /////, T23, 'R/RZERC , TAB, 'TIME TC CAVITY CC 
LLAPSE' ) 

lA FORMAT (/,T15 ,I3,T23,F8.6,TA8,F8 .6 ) 

16 FORMAT (6A8) 

STCF 

END 
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FLKCTICN FCTiX) 

f-CT = Xv«K5/( U0-X’’-*3. C)’^*C.5 

RETURN 

EEC 

//GC .SVSIN OD * 

HGLF’.D. TIME TO CAVITY CCLLAPSE FRCM RZERO 
TO A RADIUS R (PL.CTTED AS R/RZERO) 
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//HCLI-' JCB (0738, 1023, AR22) , 'hOLM' ,TIFE = /i 
// t.NEC r-CRTCLGP,RcGI0N.GC=100K 
//FCFT.SVSIN CD 



CAVITY FROBLFN MUREuR 2 



^ ?}t >{c 5^ ?}: :{« .%‘S :^c * i: ^ *:{» >) ^4 A ^ ^ 

a 

'T 

* THIS PRCBLEM USES THE .NUMERICAL INTEGRATICIx SUB- 
^ RCUTINE QATR TO SCLVE THE NCR-LINEAR CAVITY ECUATICK * 
CF MOTICK. PRCBLEM .NUMBER 2 IS CESIGNED TC CCMPUTE 
CSCILLATCRY MQTIGN QF A GAS FILLEC HYDRAULIC RAM 
^ CAVITY. THE OUTPUT IS A PLOT OF THE CAVITY PRESSURE ” 
^RATIGVSTIf-IE. 

■it ;k 

EXTERNAL FCT,F 

CIHENSICN T(500),P(500) ,AUX(500) 

CCMMON GAMMA, ALPHA 
WRITE (6,19) 

THE INITIAL CONDITIONS FOR PRESSURE AND FFCFERTIES 
CF THE GAS ARE DEFINED. 

GAMMA = A .0/3. C 
FRATI0=6.667 

•TIME' IS THE CUMULATIVE TIME CCLNTER. 

TIME=O.C 

ITMAX=5C 

N IS THE SOLUTION ARRAY COUNTER. 

N= 1 

THE FIRST ROOT OF THE CAVITY EQUATICN OF MCTICN IS 
FCLNO USING SLBRGGTINE ZFALSE. 

1 >L=C.30 
XR=C. SO 

ALFhA = PRATIC/(GAM!MA-l .0) 

CALL ZFALSE . ( F , 0 . 00 1 , 3 , XL , XR , XA P P , I TMAX , I E R ) 

FF = FRAT IC 

P RCCT=X A PP** (3.0-GAMMA ) 

PFATIO=FRCGT*PRAT IC 

INTEGRATION PARAMETERS FOR SUEROLTINE QATF ARE 
ASSIGNED . 

XLCW=0. 999 

XL=C.999 

DC 2 I=N,5C0 

IF { I. EC. 311 ) GO TC 10 

CALL QATP(XLCVI,XU,C.0C1,40,FCT,Y,IER,AUX) 

T ( I J=Y + TIME 

P( I ) = XLCU=*^(3.0*GAMMA)*FP 
XLCVx=XLCN-C. 020 
IF (XLOW.LT.XAPP) GC TC 3 

2 CCNTlNUE 

3 TIMF=T(I) 

N= I 

WRITE (6,20) I ,XAPP,T (I ) ,P( I) ,PP 
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c 

c 

c 

c 



c 

c 



c 

c 

c 

c 



c 

c 



c 

c 

c 

c 



c 

c 



XL=2.20 

Tf-E SECCf\0 ROOT OF THE CAVITY ECCATION CF RCTICN IS 
FCLf\D USIKG SLBRGLTINE ZFALSE. 

XR-2.^iO 

ALFhA=PRAT 10/ (GARRA-1 .0 ) 

CALL ZFALSE ( F , 0 . 00 1 , 3 , XL , X R , XA P P , 1 TMAX v I E R ) 

PF = PRAT IC 

FRCCT = XAPP** (3,0*GAPPA ) 

PFATI0=PRCCT-PRATIC 

IKTEGRATICfv PARAMETERS FOR SLEROUTINE QATP ARE 
ASSIGNEE. 

XLCW= 1. CCl 

XL=i.OOI 

CC 6 I=K,500 

IF (I. EG. 311) GO TC 10 

CALL QATRI XLOWrXU ,C .001 FCT, Y , lER , AUX ) 

T ( I ) = Y + TIME 

P( I ) = XU»!'*(3..0*GAMMAI«PP 
XC = XL+C .020 

IF (XU.GT.XAPP) GC TO 7 
e CONTINUE 
7 lIf'E=T{I) 

N=I 

NRITE (6,20) I,XAPP,T (I ) ,P( I ) tPP 
GC TO 1 

IC NFITE (6,13) 

CC 11 M=l,310 

NRITE (t,l^+) F,T(M),P(M) 

11 CCNTINUE 

SCLUTICN ARRAY IS PREPARED FOR PLOTTING USING SUE- 
RCCTINE CRAN. 

REAL’i'8 TITLE(12) 

PEAC (5,15) TITLE 
REAL LAELE/' '/ 

GALL DRAW ( 300 , T , P , 0 , 0 , LA B L E , T I TL E , C . 0 , 0 . C , C , C , 0 , C , 6 , 8 
,0,LAST) 

12 FORMAT («!• ,./////, T25, 'NON-DIMENSIONAL T I N L ' , T 50 , ' C A V I 

TY PRESSURE/PRESSURE AT INFINITY') 

FORMAT (/,T18,I3,T27,FIC.4,T53,FiC.A) 

15 FORMAT (6A8> 

IS FORMAT ( ' 1 ' ,/// , T22 , ' I ' ,T31 , 'XAPP ' ,T43, ‘ T IM E ' ,T55 , 'P ( I 
> ‘ ,T68, ' PP* ) 

2C FORMAT ( / , 20X , I 3 , 4 ( 2X , F 10 . 4 ) ) 

STCP 

ENC 



FLNCTIGN F(X) 

COMMON GAMMA, ALPHA 

F = -{ ALPHA'i=X*=<-' (3.0*-=GAMMA) ) + {(ALPHA+1.0)^-X**3.C)-l.C 

RETURN 

E-RD 



FLNGTICN FCT(X> 

GCMM.ON GAMMA,ALPHA 

FGT= (X^*(-2.0 ) )^( - ( ALPHA^X** ( 3.0 -GAMMA) )-^((ALPHA+i.G)* 
X«=^3.G)-i.0)**(-0.5> 

RETURN 

ENC 
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c 



c 



c: 



c 



1C 



c 



1 £ 



C 

c 

c 



2C 



2£ 

C 

3C 

C 

Q C 

90CC 

9CC£ 



, C 



SIERCUTINE ZF;>LSE ( F , E PS » NS I G , X L , >■ R , X AP F , IT R /i X , I E R ) 

lER = 0 
IC = 0 

EFSP = lO.^^^C-rNSIG ) 

FXL = F (XL ) 

FXR = F(XR) 

IF (FXL=<=FXR) 15, 1C, 5 

TERf^lNAL ERRCR 

lER = 129 
GC TO 4C 

FXL CF FXR = G 

XAPP = XR 

IF (FXL .EC. 0.0) XAPP = XL 
GC TO AC 

CCPPLTE APPRCXIPATE RCCT 
XAPP = XL + FXL^-(XR-XL)/ (FXL-FXR) 

FXAPP = F(XAPP) 

IF (ABS(FXAPP) ,GT, EPS) GO TO 2C 
GC TO AC 

DETERMINE ir,FETFER APFRCXI 
LI ES BETWEEN XAPP ANC XL 
ANC XP 

IF (FXAPP^FXL .GT. C.OI GO TO 25 

XR = XAPP 

FXR = FXAPP 

GC TO 3C 

XL = XAPP 

FXL = FXAPP 

DETERPINE IF (XR-XL .GT. 
IF (XR^XL .GT. EPSP^ABS (XR) ) GO TC 35 
GC TO AC 

CONTINUE FOR ITPAX ITERAT 

IC = IC+1 

IF ( IC. LE.50) GO TC 15 
lER = 130 
ITPAX = IC 

IF (lER .NE. 0) GC TO 9000 

GC TO 9005 

CCNTINUE 

CALL UERTST ( I ER , ' ZFA L S E • ) 

RETURN 

ENC 



SUBROUTINE UERTST ( I ER , NAME ) 



D IFENSICN 
INTEGER^2 
INTEGER 
EGUIVALENCE 



ITYP(5,A) ,I3IT(A) 



NAME (3 ) 

WARN , WARE ,TERP 
( I ei T( 1 ) ,WARN ) 



PRINTR 

(IBIT(2 ) ,WARF) 







C/iTA 


ITYP 




/'WARN' 


, ‘ ING ‘ , • • , ' 


1 






* 






* WARN ' 


, ‘ ING { ' , ' WITH' , ' 


FIX' 












'TERM' 


, • INAL ‘ , ' ' , ' 


1 












1 KCN-' 


, 'DEF I ' , ' NED ' , ' 


1 








leiT 




/ 32 


,6A 


, 128, C/ 








CATA 


PRINTR 


/ 6/ 










I £P2=IER 
IF (IER2 


.GE. 


WARN ) 


GO 


TO 


5 




c 




IEP1=A 
GC TO 2C 










NON-CEF INED 






5 


IF (IER2 


.LT . 


TERM ) 


GO 


TO 


10 




c 




I ER 1=3 
GC TO 2C 










TERMINAL 






1C 


IF (IER2 


.LT. 


WARE ) 


GC 


TO 


15 




c 




IER1=2 










WARNING! WITH 


FIX) 



( IE IT 

I 
I 
I 
I 
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GC TO 2C 

C WAPMKG 

15 IERi=l 

C EXTRACT 'N* 

2C IER2=IER2-IBIT( I5R1) 

C PRIM ERROR PESSAGE 

WRITE (PRINTRr25) ( I T Y P ( I , I ER 1 ) , I = 1 , 5 J , N A P E , I E R2 
25 FCRPATI' I M S L(UERTST) « , 5A4 , X , 3A2 , AX , I 2 ) 

RETURN 
ENC 

//GC.SYSIN DO 

HCLP,D. FYORAULIC RAM CAVITY 
PRESSURE OSCILLATION 
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